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ABSTRACT 


In  this  paper  axial  vibrations  of  whirling  bars  are 
studied  using  undefornned  coordinates.  It  is  found 
that,  in  general,  the  effect  of  whirling  is  to  lower 
the  natural  frequencies  of  the  bars.  For  the  static 
case  an  interesting  result,  which  has  not  been  pre¬ 
viously  reported,  is  obtained  as  a  special  case  of  the 
dynamic  problem.  When  the  angular  velocity  of  the 
bar  approaches  certain  critical  values,  "static  res¬ 
onances"  occur  and  the  axial  displacements  every¬ 
where  in  the  bar  tend  to  become  unboundedly  large. 
These  resonances  are  explained  physically.  The 
same  problem  is  also  worked  out  using  final  or  de¬ 
formed  coordinates  and  it  is  shown  that  the  static 
resonances  cannot  come  out  of  such  an  analysis. 
Moreover,  the  rotation  does  not  have  any  effect  on 
the  natural  frequencies  when  deformed  coordinates 
are  used.  The  results  of  the  analysis  using  unde-- 
formed  coordinates  appear  to  be  more  compatible 
with  the  physics  of  the  problem. 


I.  INTRODUCTION 


The  study  of  whirling  bars  is  of  Interest  because  they  find  application 
as  machine  elements  such  as  turbines,  propellers  and  helicopter  blades. 
More  recently,  they  have  been  employed  as  antenna  elements  of  artificial 
satellites  used  for  communications  purposes.  In  this  paper  axial  vibrations 
of  whirling  bars  are  studied  using  undeformed  coordinates.  It  is  found  that, 
in  general,  the  effect  of  whirling  is  to  lower  the  natural  frequencies  of  the 
bars.  The  advantage  of  using  undeformed  or  initial  coordinates  is  that  they 
yield,  for  the  static  case,  an  interesting  result  as  a  special  case  of  the 
dynamic  problem.  This  result  for  the  static  case  has  not  been  previously 
reported.  When  the  angular  velocity  of  the  bar,  u  ,  approaches  certain 
critical  values,  static  resonances  occur  and  the  axial  displacements 
everywhere  in  the  bar  tend  to  become  unboundedly  large.  These  reso¬ 
nances  are  explained  physically.  The  same  problem  is  also  resolved  using 
final  or  deformed  coordinates.  Such  a  resonance  can  never  result  from  an 
analysis  using  deformed  coordinates  since,  for  the  axial  displacement,  u, 
to  get  arbitrarily  large,  the  deformed  coordinate  must  become  arbitrarily 
large.  Thus,  the  only  u  that  can  become  arbitrarily  large  is  the  one  for 
which  the  deformed  coordinate  becomes  large,  and  the  resonant  character 
of  the  solution  is  lost.  As  the  angular  velocity  of  the  bar  approaches  zero, 
the  frequency  equations  reduce  to  those  of  bars  that  vibrate  axially  but  are 
otherwise  at  rest.  A  comparison  of  the  static  solutions  obtained  from  the 
use  of  deformed  and  undeformed  coordinates  is  made.  Moreover,  the  use 
of  undeformed  coordinates  does  not  even  show  the  effect  of  rotation  on  the 
natural  frequencies. 

II.  DIFFERENTIAL  EQUATIONS  OF  MOTION 

Let  (X,  Y)  be  a  reference  system  fixed  in  space  and  (x,y)  a  coordinate 
system  attached  to  the  undeformed  whirling  bar  with  its  origin  coincident 


) 
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with  that  of  the  (X.  Y)  eyetem.  The  (x,  y)  coordinate  eyetem  will  then  have 
an  angular  velocity  w  (for  a  fixed  x)  which  will  be  called  the  angular  velocity 
of  the  bar.  If  the  dieplacement  vector  is  defined  ae  the  difference  between 
the  radius  vector  to  the  particle  at  time,  t,  denoted  by  R,  and  the  radius 
vector  to  the  particle  at  some  initial  time,  t^,  denoted  by  R^,  then  the 
displacement  vector,  V,  is  given  by  R-R^  and  the  acceleration  vector  by 

a  =  - y  ([1],  p.  72).  Now  V  as  it  stands  could  be  quite  large,  and  con- 

ar 


vective  terms  must  be  used  in  differentiating  it  with  respect  to  time,  since 
small  deformation  theory  cannot  be  applied  ([2],  p.  85).  This  would  result 
in  a  nonlinear  theory,  which  is  inconvenient.  To  eliminate  the  use  of  con> 
vective  terms,  it  is  assumed  that  the  particle  does  not  deform  greatly  from 
the  xmdeformed  rotating  position  (the  x,  y  system  rotates  with  the  angular 
velocity  u).  Define  a  displacement  v  such  that  V  =  r  -f  v,  where  r  is 
the  position  vector  to  the  rotating  point  located  at  x.  The  acceleration 


will  then  be  given  by 


at^ 


unit  vector  e^  in  the  x  direction  and 
angular  velocity  vector  w  is  given  by 


To  carry  out  this  process,  set  in  a 
e  in  the  y  direction.  Then  the 

y  ' 


fa)  s  ta>  e. 


(1) 


and 


r  =  (x  +  u)  e^  ,  (2) 

where  u  is  the  magnitude  of  the  component  of  v  in  the  x  direction. 

The  other  components  of  v  in  the  direction  of  the  tangent  in  the  plane 
of  rotation,  and  in  a  direction  transverse  to  the  plane  of  rotation,  are  ig¬ 
nored  since  only  the  axial  vibrations  of  whirling  bars  are  considered  in  this 
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paper.  One  of  this  paper'*  object*  i*  to  bring  out  the  advantage  of  u*ing 
undeformed  or  initial  coordinate*  in  thi*  problem.  Carrying  out  the  differen¬ 
tiation  indicated  earlier,  the  equation  of  motion,  employing  the  u*ual  a**ump- 
tion*  for  the  one-dimen*ional  elaetic  wave  propagation  in  bar*,  i* 

00-  0^u  2,  .  .  ... 

irr  -  P  ~~T  " 

*  0t 

where  <r  i*  the  normal  *tre**  in  the  x  direction  and  p  i*  the  deneity  of  the 
material.  To  obtain  equilibrium  equation*,  deformed  coordinate*  *hould  be 
u*ed.  However,  in  the  limit  of  email  etrain  theory  (which  aeeumption  i* 
made  here),  the  equilibrium  equation*  and  the  *tre**e*  become  identical 
( [1]  ,  Chap.  5)  for  both  *y*tem*  of  coordinate*.  The  *tre**-*train  relation 
is 


(T  = 


Introducing  (4)  into  (3)  yield* 


where 


2  e^u 


et‘ 


2  2 
or  u  u 


2  2 


or  u  X 


2 


a 


t 


(4) 


(5) 


(6) 


and  E  is  the  modulus  of  elasticity. 
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III.  BOUNDARY  CONDITIONS 

Two  sets  of  boundary  conditions  will  be  considered.  The  first  set  nnay 
pertain  to  the  vibrations  of  a  helicopter  blade  or  a  communications  satellite 
antenna,  whereas  the  second  set  pertains  to  the  web  of  a  flywheel  with  a 
rigid  flange  on  the  outside. 

(a)  The  boundary  conditions  for  the  first  case  are: 

u(0.t)  =  0  ,  (7a) 


= 


0 


(7b) 


which  correspond  to  the  case  of  a  fixed-free  bar. 

(b)  The  boundary  conditions  for  the  second  case  are: 


u(0,t)  =  0  , 

(8a) 

u(i,  t)  =  0  , 

(8b) 

which  correspond  to  the  case  of  a  fixed-fixed  bar. 

IV.  SOLUTIONS  OF  BOUNDARY  VALUE  PROBLEMS 
For  the  solution  of  the  boundary  value  problem,  split  u  such  that 

u(x,  t)  =  Uj(x)  +  ,  (9) 


where  u^  and  U2  satisfy  the  following 


2  2 
a  b>  u, 


2  2 


a  u  X 


(10) 
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and 


2  2 
or  u)  u. 


(11) 


The  general  solution  of  (10)  gives  the  solution  for  the  static  problem  and 
may  at  once  be  written  as 


u^  =  a^  sinorfauc-i-  a2Cosor(t)x  -  x 
A  solution  of  (11)  may  be  taken  as 

U2(x.t)  =  U{x)eK^ 


(12) 


(13) 


where  is  the  frequency  to  be  determined.  Substitution  of  (13)  into  (11) 
results  in 


^  +  ar^(«^+p2)  U  =  0 
dx^  ” 


The  general  solution  of  (14)  is 


(14) 


U  = 


n  T”  n — r 

a^sinor/w+p^  x  +  s^coso/w  +p^ 


(15) 


Case  (a):  fixed-free  bar 

Since  the  boundary  conditions  (7)  have  to  be  satisfied  for  all  values  of 
time,  it  follows  that  u^  and  U2  each  have  to  satisfy  (7).  It  may  be 
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remarked  that  the  static  solution  for  as  given  by  (12),  differs  from  the 
one  given  by  Biezeno  and  Grammel  ([3],  p.  124)  in  that  final  coordinates 
are  used  in  [3]  ,  whereas  the  use  of  initial  coordinates  in  the  present  paper 
has  resulted  in  a  different  differential  equation.  Application  of  (7a)  and 
(7b)  results  in 


a 


2  "" 


0 


(16) 


and 


a 


1 


1 

oru  cos  oKtil 


(17) 


provided  that  ui  does  not  correspond  to  a  root  of 

cosorui  =  0  (18) 


For  the  present  it  will  be  assumed  that  u  does  not  correspond  to  a  root  of 
(18).  The  case  when  u  corresponds  to  a  root  of  (18)  will  be  considered 
separately.  Hence  the  static  solution  is 


u 


1 


singciBc 
erucos  aul 


X 


(19) 


For  the  case  when  <•>  corresponds  to  one  of  the  roots  of  (18),  a  solution  to 
(10),  satisfying  the  boundary  conditions  (7a)  and  (7b),  does  not  exist.  When 
one  encounters  such  nonexistence  of  a  solution  to  a  physical  problem  it 
must  have  significance.  The  physical  problem  must  have  some  sort  of  a 
solution  even  though  it  only  tends  toward  infinity  as  w  approaches  w,  one 
of  the  critical  values  [roots  of  (18))  .  As  w  approaches  the  constant  a^ 
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given  by  (17)  becomes  very  large.  Such  values  o£  are  given  by 

"i  =  (^)f^  - •  <“> 

Thus,  the  lowest  angular  velocity  at  which  the  static  resonance  occurs  in  a 
£ixed>free  bar  is  the  velocity  at  which  the  peripheral  speed  o£  the  bar  is 
equal  to  times  the  speed  o£  sound  in  the  material  o£  the  bar.  At  such 
speeds  the  solution  £or  u^  becomes  large  £or  all  values  o£  x  between  aero 
andi.  Hence,  one  can  say  that,  £or  any  value  o£  x  between  aero  andi  , 

Uj  becomes  unboundedly  large  as  u  approaches  This  is  evidently  some 
sort  o£  a  resonance  phenomenon  and  is  directly  a  result  o£  the  use  o£  unde¬ 
formed  coordinates  in  the  present  analysis.  At  first  sight  it  would  appear 
that  u^  cannot  become  arbitrarily  large,  since  there  is  no  steady  energy 
input  to  make  it  do  so.  However,  this  reasoning  is  fallacious.  Consider 
the  bar  as  being  undeformed  but  having  an  angular  velocity  w  at  some  in¬ 
stant  of  time  t^.  Suppose  that  at  this  time  the  bar  is  allowed  to  deform. 

The  axial  displacement  for  the  steady  state  or  static  case,  Up  will  im¬ 
mediately  begin  to  increase,  hopefully  to  attain  some  steady  value.  As  u^ 
increases,  the  moment  of  inertia  will  increase,  since  the  elemental  mass 
will  now  be  located  at  a  distance  x+u^  instead  of  at  x.  The  angular  mo¬ 
mentum  of  the  bar  will  then  increase,  since  u  is  assumed  to  be  constant, 
and  energy  must  be  added  to  the  system.  It  is  not  unreasonable  to  suppose 
that  there  will  be  some  value  of  w  such  that  this  process  does  not  lead  to 
a  steady  bounded  value  of  u^,  and  the  roots  of  (18)  are  just  such  values  of 
w.  Actually,  only  the  lowest  value  of  is  of  consequence  since  it  is  un¬ 
likely  that  equilibrium  physically  can  be  attained  for  w  greater  than  the 
lowest  value  of 
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For  u^i  the  application  of  boundary  conditions  (7a)  and  (7b)  to  (15) 
results  in 


a 


4  " 


0 


(21) 


and 


COS 


aj} 


27T  * 


For  a  nontrivial  solution  it  is  required  that 


(22) 


aj^ 


2^  2  . 
+Pn^  = 


(2n+l)  n=0.1.2, 


(23) 


From  (23)  the  frequency  equation  can  be  rewritten  as 


n 


n=0, 1,2, 


(24) 


If  u  is  allowed  to  approach  zero,  (24)  reduces  to  the  frequency  equation  of 
a  fixed-free  bar  that  is  vibrating  axially,  but  is  otherwise  at  rest.  Equation 
(24)  shows  that,  in  general,  the  effect  of  whirling  is  to  lower  the  natural  fre¬ 
quencies.  It  is  of  interest  to  consider  two  special  cases  arising  from(24) 
as  w  increases.  First,  it  is  possible  for  w  to  be  large  enough  so  that  the 
right  member  of  (24)  vanishes  for  a  certain  value  of  n.  It  can  easily  be 
shown  by  referring  to(15)that  this  corresponds  to  a  static  instability  in  the 
particular  n^^  mode.  It  should  be  remarked  that  this  static  instability  is  the 
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same  as  the  one  which  has  already  been  mentioned  for  u^.  Second,  consider 
the  case  when  u  is  allowed  to  be  large  enough  so  that  the  right  member  of 
(24)  becomes  negative  for  values  of  n  <  n^  and  is  positive  forn>n^.  Then 
for  n>n  ,  vibrations  exist  whereas,  for  n<n  ,  p  becomes  imaginary  and 

C  C  II  ^  * 

these  modes  could  get  damped  out  or  grow  asymptotically  with  time  depend¬ 
ing  on  the  sign.  However,  it  is  not  meaningful  to  consider  the  vibrations  of 
the  bar  in  this  domain,  since  it  is  unlikely  that  static  equilibrium  can  be 
attained  for  such  values  of  w.. 

Case  (b):  fixed-fixed  bar 

Application  of  boundary  conditions  (8)  to  (12)  gives 

a2  =  0  ,  (25) 


and 


i 

~  sin  aui  ’ 

provided  that  w  does  not  correspond  to  a  root  of 

sin  owi  =  0 


(26) 


(27) 


Hence, 


u 


1  ~  slrTowJ 


sinouK  -  X 


When  u  corresponds  to  a  root  of  (27);  i.  e.  , 


iir  IE  •  1  1  ■j 

"i -tJj  • 


(28) 


(29) 
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static  resonances  occur,  as  has  already  been  explained  in  the  case  of  a 
fixed-free  bar.  The  lowest  angular  velocity  at  which  the  static  resonance 
occurs  in  a  fixed-fixed  bar  is  the  velocity  at  which  the  peripheral  speed  of 

t 

the  bar  is  equal  to  ir  times  the  speed  of  sound  in  the  material  of  the  bar. 
Only  the  lowest  value  of  uj.  is  of  interest  as  has  been  remarked  in  the  case 
of  a  fixed-free  bar. 

Similarly,  when  applied  to  u^,  the  boundary  conditions  (8a)  and  (8b) 

yield 


^4 


=  0 


(30) 


and 


a^sin 


/ITT 

in  +p^ 


i  =  0 


(31) 


For  a  nontrivial  solution  it  is  required  that 


rT~T 

V"  *”1. 


/  -  nir  ,  n=l ,  2,  3, . . . 


(32) 


from  which  the  frequencies  for  the  fixed-fixed  bar  are 


2  _  n^ir^  E  2  ,  _  _ 

p  —  •  ca>  »  nsj 

"  ^2  p 


(33) 


If  u  approaches  zero  in  (33),  the  resulting  equation  for  p^  is  the  same  as 
the  one  for  a  fixed-fixed  bar  that  is  vibrating  axially  but  is  otherwise  at 
rest.  Remarks  pertinent  to  the  ones  for  the  fixed-free  bar  apply  also  here, 
when  (Ji  becomes  successively  large. 
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V.  SOLUTION  USING  DEFORMED  COORDINATES 


If  now  deformed  or  final  coordinates  are  used  in  obtaining  the  inertial 
forces,  the  equation  of  motion  (5)  becomes 


2 


a 


2 

dr  u  X 


(34) 


Split  u  such  that 


u(x,t)  =  Uj(x)  +  U2(x,t) 


(35) 


where  and  U2  satisfy  the  following 


+ 


2  2 


or  w  X 


0 


and 


(36) 


(37) 


The  general  solution  of  (36)  is  the  solution  for  the  static  problem  and  may 
be  written  as 


u. 


2  2  3 
a  w  X 


(38) 
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Equation  (38)  is  identical  to  the  one  that  can  be  deduced  from  [3],  For  the 
fixed-free  bar  the  static  solution,  after  evaluating  the  constants  of  integra¬ 
tion  using  (7a)  and  (7b),  is 


“l  =  —I—  (  '  •  T  )  ■  (”> 

Equation  (39)  shows  that  the  resonant  character  of  the  static  solution  is 
indeed  lost  when  deformed  coordinates  are  used.  Furthermore,  it  follows 
from  (37)  that  the  frequency  equation  for  a  whirling  bar  that  is  vibrating 
axially  would  be  independent  of  the  angular  velocity,  u  ,  with  which  the  bar 
is  whirled.  A  steady  input  of  energy  can  be  made  into  the  system  free  of 
any  dissipation  without  an  instability  ever  occurring,  a  result  which  does 
not  appear  to  be  realizable. 

VI.  NUMERICAL  EXAMPLE 

To  illustrate  the  results  of  the  analysis,  numerical  examples  were 
worked  out  for  a  steel  bar  100  inches  in  length.  The  static  displacements 
for  a  fixed-free  bar,  as  given  by  deformed  and  undeformed  coordinates, 
are  given  in  Fig.  1.  The  ratio  of  the  displacement  at  the  tip  of  the  fixed- 
free  bar  obtained  from  the  use  of  undeformed  and  deformed  coordinates  is 
plotted  as  a  function  of  u  in  Fig.  2.  The  example  shows  that  for  small 
values  of  w  the  two  results  do  not  differ  appreciably,  but  the  results  of  the 
analysis  using  undeformed  coordinates  indicate  a  new  limit  to  the  speed  at 
which  the  bars  can  be  whirled  even  with  an  exotic  material  having  a  very 
high  yield  point.  Figures  3  and  4  [obtained  from  (24)  and  (33.)]  show  the 
effect  of  whirling  on  the  natural  frequencies  of  fixed-free  and  fixed-fixed 
bars,  respectively.  The  intersection  of  the  curve  and  the  abscissa 
indicates  the  occurrence  of  instability. 
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UNOEFORMEO  /  u,  DEFORMED 


40 


L*  too  IN. 
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Fig.  3  -  Plots  of  vs  w  for  a  Fixed-Free  Bar 
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Fig.  4  -  Plots  of  vs  h)  for  a  Fixed-Fixed  Bar 
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vn.  CONCLUSION 


The  effect  of  whirling  on  the  axial  vibrations  of  a  bar  has  been  studied 
in  this  paper.  Certain  advantages  of  the  use  of  undefornned  coordinates  in 
this  problem  are  indicated  and  the  results  of  such  an  analysis  appear  to  be 
more  compatible  with  the  ]^ysics  of  the  problem.  When  the  angular  velocity 
of  the  bar  reaches  certain  critical  values,  static  resonances  occur.  These 
resonances  point  to  a  new  limit  to  which  the  bars  can  be  whirled  even  with 
the  use  of  an  exotic  material  having  a  very  high  yield  point. 
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